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Ultrametric structure of the energy cascade process in a dy- 
namical model of turbulence is studied. The tree model we use 
can be viewed as an approximated one-dimensional trunca- 
tion of the wavelets resolved Navier-Stokes dynamics. Vary- 
ing the tree connectiveness, the appearance of a scaling tran- 
sition in the two-points moments of energy dissipation is de- 
tected, in agreement with experimental turbulent data. 



Spatio-temporal intermittency is the most intriguing 
aspect of a fully developed three-dimensional turbulent 
flow. 

Experiments 0] show that the energy dissipation de- 
flnes a multifractal measure on the fluid volume. The 
multifractal measure is characterized by the scaling prop- 
erties of the coarse-grained energy dissipation on a box 
at scale r, Sr, namely: < {erix)y >~ r'^'-P\ where < • > 
means averaging over all boxes of size r and centered in 
X in which the volume can be partitioned. The measured 
t{p) exponents show a clear intermittent behaviour, i.e., 
a non-linear dependency on the order of the moment p. 

The simplest way to explain phenomenologically the 
presence of intermittent deviations consists in describing 
the energy transfer mechanism in terms of fragmenta- 
tion stochastic processes. In these models (see for a 
recent proposal), one introduces a set of eddies leaving 
on a dyadic structure and connected through a random 
multiplicative process. 

Let us remark that all stochastic fragmentation models 
so far proposed lack any direct linking with the original 
Navier-Stokes (NS) equations. Dynamical deterministic 
models on hierarchical structures are therefore invoked 
for improving our understanding of the energy transfer 
mechanism. 

In this paper, we study dynamical models which fill the 
gap between purely stochastic fragmentation models and 
the original NS dynamics. In particular, we consider a 
dynamical model on one-space and one-time dimensions 
1^. One can look at this model as an approximation of 
the original NS equations in a wavelets basis |^ . 

In order to specify the model one has to select the set 
of interactions connecting eddies at different scales and at 
different spatial locations. By changing the interactions 
set, one changes the scale-organization of energy struc- 



tures: in order to study it, new tools are required 
which characterize intermittency more completely than 
the multifractal spectrum r(p) alone. An obvious gener- 
alization of the single-point statistics is to inquire about 
the scaling of two-point moments. In practice, introduc- 
ing the mixed moments; 



(1) 



one can study correlations between different scales by 
changing r and / and/or correlations between different 
points in the fluid volume by changing s. 

Different interactions among nodes of the hierarchical 
dynamical structure can lead to very different prediction 
for the scaling behaviour of (Q). An ultrametric organi- 
zation of the tree can be detected by looking for a phase 
transition in the set of scaling exponents characterizing 
correlations (^. In the analysis of the two-points ob- 
servables (0) performed on experimental turbulent data 
gave a first support for an ultrametric organization of the 
main triadic interactions in Navier-Stokes eqs. 

In the following, we are going to analyse the same kind 
of observables measured on a direct numerical integra- 
tion of a dyadic-tree model for turbulent energy cascade. 
In particular, by changing the set of interacting triads 
we want to disentangle the basic symmetries behind the 
transition observed in real turbulent data. 

Let us turn to a brief review of the model (for a com- 
prehensive description see |^). The tree model can be 
viewed as an extension of shell models, which can be seen 
as a severe truncation of the NS equations (see ^ for a 
general introduction). The most popular shell model is 
the Gledzer-Ohkitani-Yamada (GOY) model ( §- Q). 
Recently, a new class of shell models based upon the he- 
lical decomposition of NS equations [|l^ has been sug- 
gested and studied ]l5| , p^ . In these models, one or 
few complex variables w„ represent an entire shell of wave 
numbers k such that kn < k < kn+i, with kn =2". Shell 
models can be thought of as field problems in zero spatial 
dimension (d=0). In order to include also some real space 
dynamics we need to transform the chain-model into a 
tree-model with d ^ 1. This is achieved by letting grow 
the number of degrees of freedom with the shell index n 
as 2". The tree model can be regarded as describing the 
evolution of the coefficients of an orthonormal wavelets 
expansion of a one-dimensional projection of the veloc- 
ity field. In the tree model, we use the notation ^ to 
indicate a complex variable having positive or negative 



1 



defined helicity and living on scale kn and spatial posi- 
tion labeled by the index j . For a given shell n, the index 
j can vary from 1 to 2"^^. 

We report here the structure of the tree model dynam- 
ical equations (for more details, see ^): 

-iyklul^+d„^noF+. (2) 

Here, n = 1,...,N, where N is the total number of 
shells, v is the viscosity, the external forcing act- 
ing on a the large-scale shell hq = 1 and ani,n2,ji,j2 
parameters, which are determined by imposing conserva- 
tion of energy and helicity in the inviscid and unforced 
limit. The si and S2 indices are the helicity signs (±) 
of interacting modes. The same equations hold, with all 
helicities reversed, for u~j. 

In restricting the possible choices of the nonlinear 
terms, we can phenomenologically require a certain de- 
gree of locality for interactions among variables at differ- 
ent scales and at different spatial positions. Regarding 
scale numbers ni and n2, in our equations (H) each vari- 
able u'^ J is allowed to interact only with nearest and 
next-to-nearest levels: indeed, ni and n2 can vary only 
from (n — 2) to (n -|- 2). Regarding space numbers ji 
and j2, we define two different models having different 
topological structure of the dynamical interactions. The 
first model, hereafter called model A, is a model which 
has an ultrametric structure. That is, each eddy is al- 
lowed to interact only with bigger eddies which spatially 
contain it and with smaller eddies spatially contained 
in it. This model is the natural dynamical representa- 
tive of all stochastic fragmentation models which phe- 
nomenologically reproduce single-point intermittent ex- 
ponents. The second model, hereafter called model B, 
has an enlarged interactions set containing also horizon- 
tal couplings, which allow eddies covering different spa- 
tial regions to interact each other. In figures 1 and 2, we 
pictorially show the set of interactions defining models A 
and B. 

The single-point statistical properties of the tree 
model have been studied in ||]: in both cases A and 
B, the system turned out to have an intermittent energy 
transfer qualitatively similar to what one can find in the 
original NS eqs. The tree-like structure imposed on the 
velocity fluctuations does not necessary implies that the 
energy dissipation can be described in terms of fragmen- 
tation processes. In order to test the scale-organization of 
the energy structures, ultrametric-sensitive observables 
should be studied. 

In order to compare our system with previous findings 
we shall detect possible ultrametric structure in the 
energy cascade of the tree model looking at two-point 
statistical quantities. All parameters settings and nu- 
merical methods are as in 0. In particular, we consider 



a total number of levels N = 16: the total number of 
sites forming the tree is then Nt = 2^ — 1 = 65535, each 
one described in terms of two complex variables. Numer- 
ical simulations needed state-of-the-art multi-processor 
computers. 

The fields we focus on are the coarse-grained energy 
dissipation densities, here denoted as e„(j), obtained as 
averages over spatial regions Aj(n) of length 2"". We 
consider the mixed moments (|l|) with r — I and p — —q, 
which in our notation become: 

<el{j)e-\j+s)>^Cn{q.s) (3) 

The behaviour of this quantity for intermittent ultra- 
metric measures resulting from random multiplicative 
processes has been already analyzed in the framework of 
the two-point multifractal formalism In this case, 

the average can be properly decomposed and a general 
result can be obtained for its dependence on the spatial 
distance s between the two points: 

This expression implies that for some moment q, a 
sharp transition occurs in the derivative of the scaling 
exponent ^(q). This scaling transition is the analog of 
a phase transition in the thermodynamic interpretation 
of multifractals ||l8|. The behaviour of (||) is dominated 
by pairs of points at which the dissipation is very large 
at one point and very low at the other. This constitutes 
the subset of points that are likely to be independent 
from each other and lie on the boundary of their bigger 
predecessors. Indeed, it must be recalled that in a ultra- 
metric structure nearby (in space) eddies could lie on the 
boundaries of much bigger ones, then having an effective 
large ultrametric distance. 

The two spatial scales of interest are the coarse- 
graining scale /„ = 2~" and the offset scale Is — Z„s: 
they should be such that rj <^ In <Si h Ay, where rj 
and At are the Kolmogorov and integral scale, respec- 
tively. For this reason, we fixed n = 11, in order to 
consider the largest inertial scale in our tree structure, 
and we let s to assume the exponentially spaced values 
s — 2™, with m ~ 1,2,. ..,71 — 2. In order to test the 
presence of a scaling transition, the mixed moments (|^) 
have been computed for 0.5 < g < 4. 

Figure 3 shows the mixed moments (^ as a function 
of s in log-log coordinates, for increasing values of q and 
for the two versions A and B of the tree model. The ^{q) 
exponents have been calculated by linear fit in the inertial 
range region: they are reported in figure 4, where are 
compared with the curve {—T{q) — T{—q)), obtained using 
the single-point moments exponents and corresponding 
to the predicted form of $(g) if the scaling transition were 
absent. In the case A, the data support a sharp transition 
in the derivative of ^{q) at q 1.5, with a much slower 
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variation of ^{q) for q > 1.5. This transition is absent in 
case B. 

We thus conclude that version A gives support for 
a scaling transition in the mixed moments of coarse- 
grained dissipation. This result agrees with the experi- 
mental behaviour found in ||^ (see fig. 17 of this reference) 
using data measured in a turbulent wake. The physical 
picture implied by this scaling transition is that of uncor- 
related small eddies that come close together even sharing 
no common history during the energy cascade. 

Let us summarize our results. 

A dynamical model in one spatial dimension orig- 
inating from a wavelets-like decomposition of a one- 
dimensional cut of a turbulent velocity field has been 
studied. We found that a scaling transition appears as 
soon as the tree has a pure ultrametric dynamical struc- 
ture. 

The fact that decreasing the number of triad interac- 
tions one can reproduce the real data scaling transition 
observed in |^ may seem in contrast with the observa- 
tion that in the original Navier-Stokes equations all pos- 
sible interactions are switched on. This contradiction is 
only apparent: divergenceless character of the original 
NS field, added to complex phase-coherence effects, can 
very easily introduce different dynamical weights in the 
possible triad interactions leading to a situation where 
only a few of them govern the global dynamical evolution. 
For example, Grossmann and coworkers showed p9| , p0[ , 
by performing suitable truncation of NS equations, that 
intermittency depends on the typical degree of locality 
in Fourier space of the survived triad interactions; very 
similar results have also been found in shell models at 
varying the inter-shell ratio A [ p6[ . 

Recently, some theoretical studies and experimental 
analysis have been done on multi-point multi-scale ve- 
locity correlation functions in turbulent flows j2|]. Our 
dynamical investigation suggests that, in the presence of 
a strong ultrametric structure, correlations among veloc- 
ity fluctuations at different scales should depend on their 
ultrametric distance rather than on the separation length 
only, as predicted in ||2l|| . 

These and similar studies performed on such kind of 
models can improve our understanding of basic mecha- 
nisms underlying turbulent cascade. For instance, it is 
important to recognize those interactions which are more 
effective in the energy transfer mechanism when con- 
structing eddy-viscosity models and in simulating small 
scale statistics by some closure approach. 
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Figures Captions 



Figure 1: Pictorial representation of interactions in 
model A: in this case the nonlinear terms in (H) 
result from the sum of parts (a) , (b) and (c) in the 
figure. 

Figure 2: Pictorial representation of interactions in 
model B: in this case the nonlinear terms in (|^) 
result from the sum of parts (a), (b), (c), (d), (e) 
and (f) in the figure. 

Figure 3: Log-log plot of mixed moments of order 
q = 0.5, 1, 1.5, 4 (from bottom to top) against 
the distance s, for model A (left)and B (right). 
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- Figure 4: The exponents ^{q) (black circles) as a 
function of q for model A (left) and B (right). For 
comparison, the values of the function {—T{q) — 
T{—q)) are also reported (white circles). 
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